Experimental demonstration of the topological surface states protected by the 

time-reversal symmetry 
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We report direct imaging of standing waves of the nontrivial surface states of topological insula- 
tor Bi2Te3 by using a low temperature scanning tunneling microscope. The interference fringes are 
caused by the scattering of the topological states off Ag impurities and step edges on the Bi2Te3(lll) 
surface. By studying the voltage-dependent standing wave patterns, we determine the energy dis- 
persion E{k), which confirms the Dirac cone structure of the topological states. We further show 
that, very different from the conventional surface states, the backscattering of the topological states 
by nonmagnetic impurities is completely suppressed. The absence of backscattering is a spectacular 
manifestation of the time-reversal symmetry, which offers a direct proof of the topological nature of 
the surface states. 

PACS numbers: 68.37.Ef, 73.20.-r, 72.10.Fk, 72.25.Dc 



The strong spin-orbital coupling in a certain class of 
materials gives rise to the novel topological insulators in 
two 1, 2] and three dimensions [1, Ij, S @, 0] in the 
absence of an external magnetic field. The topological 
states on the surfaces of three dimensional (3D) mate- 
rials have been studied recently in Bii__a;Sb2; [1, H, 0, 
BizTea and BiaSeg 0, E, [H E El, which possess 
insulating gaps in the bulk and gapless states on sur- 
faces. The surface states of a 3D topological insulator 
comprise of an odd number of massless Dirac cones and 
the crossing of two dispersion branches with opposite 
spins is fully protected by the time-reversal-symmetry 
at the Dirac points. Such spin-helical states are ex- 
pected to bring forward exotic physics, such as magnetic 
monopole |14i| and Majorana fermions J^, JJ)]. To date, 
the experimental study of topological insulators is pre- 
dominantly limited to the determination of their band 
structure by angle-resolve photoemission spectroscopy 
(ARPES) % i, M Ei E 111 ■ Distinct quantum phe- 



nomena associated with the nontrivial topological elec- 
tronic states still remain unexplored. Particularly, there 
is no direct experimental evidence for the time reversal 
symmetry that protects the topological property. Here, 
using the low temperature scanning tunneling microscopy 
(STM) and spectroscopy (STS), we report the direct ob- 
servation of quantum interference caused by scattering of 
the 2D topologically nontrivial surface states off impuri- 
ties and surface steps. Our work strongly supports the 
surface nature of the topological states, and provides a 
way to study the spinor wave function of the topological 
state. More significantly, we find that the backscattering 
of topological states by a nonmagnetic impurity is forbid- 
den. This result directly demonstrates that the surface 
states are indeed quantum-mechanically protected by the 
time reversal symmetry. 



Iig . 1201 1 result from the 2D surface states perturbed by 
surface defects. A surface state is uniquely characterized 
by a 2D Bloch wave vector k within the surface Brillouin 
zone (SBZ). During elastic scattering, a defect scatters 
the incident wave with a wave vector ki into kf = ki + q, 
with ki and kj being on the same constant-energy con- 
tour (CEC). The quantum interference between the ini- 
tial and final states results in a standing wave pattern 
whose spatial period is given by 27r/q . When the STM 



images of a standing wave are Fourier transformed 21|, 
the scattering wave vector q {hq is the momentum trans- 
fer) becomes directly visible in the reciprocal space. In 
contrast, for bulk states, there will be continuous ranges 
of wave vectors on the projected SBZ for a given energy. 
Usually, no distinct interference fringe can be produced 
by bulk states and visualized by STM. In this sense, the 
standing wave is surface-states-sensitive and particularly 
suitable for studying topological insulators. 

Our experiments were conducted in an ultra-high vac- 
uum low temperature (down to 0.4 K) STM system 
equipped with molecular beam epitaxy (MBE) for film 
growth (Unisoku). The stoichiometric Bi2Te3 film, a ro- 
bust topological insulator, was prepared on single crystal 
substrate Si(lll) by MBE. Details of sample prepara- 



The interference patterns in STM experiments [17|, |18 



tion are described elsewhere [22]. Shown in Fig. 1(a) 
is a typical STM image of the Bi2Te3 film with a thick- 
ness of ~100 nm. The atomically flat morphology of 
the film is clearly observed. The three steps seen in 
Fig. 1(a) all have the height (0.94 nm) of a quintuple 
layer. The steps are preferentially oriented along the 
three close-packing ([100], [110] and [010]) directions. 
The image with atomic resolution [Fig. 1(b)] exhibits 
the two-dimensional hexagonal lattice structure of the 
Te-terminated (111) surface of Bi2Te3. Our STM obser- 
vation further reveals a small density of clovershaped de- 
fects on the surface (see supporting material [11). Simi- 
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FIG. 1: (a) The STM topograph of the Bi2Te3(lll) film. The 
imaged area, 250 nm by 250 nm, was scanned at a sample 
bias of 3 V and tunneling current of 50 pA. (b) The atomic- 
resolution image. Tellurium atom (pink colored) spacing is 
about 4.3 A. The image was scanned at a sample bias of - 
40 mV and tunneling current of 0.1 uA. (c) dl /dV spectrum 
taken on bare Bi2Te3(lll) surface. The spectrum was mea- 
sured with setpoint V=0.3 V, 1=0.1 uA. The arrows indicate 
the bottom of conduction band and the top of valence band, 
respectively, (d) Calculated band structure of Bi2Te3(lll) 
along high-symmetry directions of SBZ (see the insert). The 
red regions indicate bulk energy bands and the purple regions 
indicate bulk energy gaps. The surface states are red lines 
around the V point. 



lar to Bi2Se3 3) 24, 25|, these structures can be assigned 
to the substitutional Bi defects at the Te sites by examin- 
ing their registration with respect to the 1 x 1-Te lattice 
in the topmost layer. 

The surface states of Bi2Te3 were investigated by STS 
and the first-principles calculations 0. The STS de- 
tects the differential tunneling conductance dl/dV [Fig. 
1(c)], which gives a measure of the local density of states 
(LDOS) of electrons at energy eV. The two shoulders in- 
dicated by arrows in Fig. 1(c) are the bottom of the 
bulk conduction band and the top of the bulk valence 
band, respectively. The Fermi level (zero bias) is within 
the energy gap, indicating that the film is an intrin- 
sic bulk insulator jS^]. The differential conductance in 
the bulk insulating gap linearly depends on the bias and 
is attributed to the gapless surface states. These fea- 
tures in STS are in good agreement with those obtained 
by the first-principles calculations (see supporting mate- 



rial [23|). According to the calculations [Fig. 1(d)], the 
topological states of Bi2Te3 form a single Dirac cone at 
the center (f point) of the SBZ [l^, 11 1, giving rise to 
a vanishing DOS in the vicinity of fc = 0. However, the 
surface states around F point overlap in energy with the 
bulk valence band. For this reason, the Dirac point is 



invisible in STS. 

On the aforementioned surface, we deposited a small 
amount (0.01 ML) of Ag atoms, which form trimmers on 
the surface, as shown in Fig. 2(a) and more clearly in 
supporting material [1^. The atomically resolved STM 
image (see supporting material [1^) reveals that the Ag 
atom in a trimmer adsorbs on the top site of a surface 
Te atom [26j. This situation is schematically shown in 
Fig. 2(b). In addition, the Fermi level shifts upwards 
in energy by 20~30 meV after Ag deposition, suggest- 
ing electron transfer from the Ag atoms to the substrate 
(see supporting material [3]). The dl/dV mapping was 
then carried out in a region containing Ag trimmers. At 
each data point, the feedback was turned off and the bias 
modulation was turned on to record dl/dV. This proce- 
dure resulted in a series of spatial mapping of LDOS at 
various bias voltages. 

Figures 2(c) to 2(g) summarize the dl/dV maps for 
bias voltages ranging from 50 mV to 400 mV from the 
area shown in Fig. 2(a). The first striking aspect of these 
images is the existence of standing wave [ij, 18, [3| in 
the vicinity of the Ag trimmers. The spatial modulation 
of LDOS by an Ag trimmer forms a hexagonal pattern, 
whose edges are perpendicular to the f — M directions in 
SBZ. This situation is more clearly resolved at large bias 
voltages [Figs. 2(f) and 2(g)]. As expected, the interfer- 
ence pat tern is anisotropic as a result of the hexagram 
CEC [10| . The spatial period of the standing wave scales 
inversely with the bias voltage and is determined by the 
momentum transfer during scattering at a given energy. 
Below 50 meV, the fringes become obscured. It results 
from a combination of two effects: (i) The wavelength in- 
creases rapidly as the bias voltage approaches the Dirac 
point, where k = 0. (ii) At low energy, more topological 
surface states with different wavelengths are involved in 
the formation of standing wave as indicated by the first- 
principles calculations [Fig. 1(d)]. The superposition of 
waves with various wavelengths smears out the interfer- 
ence fringes. With increasing bias, especially after the 
surface states in the F — M direction merge into the bulk 
conduction band at ^0.2 eV above the Dirac point ac- 
cording to calculation [Fig. 1(d)], the contribution of 
states in the f — M direction vanishes and the states in 
the vicinity of F— K direction gradually gain more weight, 
leading to more distinct interference patterns. After the 
surface states in the F — K direction merge into the bulk 
conduction band at ~0.6 eV above the Dirac point [Fig. 
1(d)], the standing waves fade out again. 

To quantify the standing waves and obtain the scat- 
tering wave vectors, we performed Fourier transforma- 
tion of the dl/dV maps into the g-space [Figs. 2(h) to 
2(1)]. One important feature in the power spectra can be 
immediately discerned by comparing the six-fold sym- 
metric pattern in the g-space with SBZ (the red hexagon 
in Fig. 2(h)): the regions with high intensity are al- 
ways oriented toward the F — M directions, while the 
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FIG. 2: Standing waves induced by Ag trimmers on Bi2Te3(lll) surface, (a) STM image (28 nm by 28 nm) of a region with 
four Ag trimmers adsorbed on Bi2Te3(lll) surface, (b) The adsorption geometry of Ag trimmer, (c) to (g) The dl /dV maps 
of the same area as (a) at various sample bias voltages. Imaging conditions: 1=0.1 nA. Each map has 128 by 128 pixels and 
took two hours to complete. The interference fringes are evident in the images. The green and red regions indicate modulations 
with high intensity, (h) to (k) The FFT power spectra of the dl /dV maps in (c) to (g). The SBZ in (h) is superimposed on 
the power spectra to indicate the directions in g-space. The resolution of FFT, which is 27r/28 nm~^, is determined by the size 
of the STM image. 



intensity in the F — K directions vanishes (see support- 
ing material [23|). Such phenomena can be understood 
by exploring possible scattering processes on the CEC 
in the reciprocal space [Fig. 3(a)]. Generally, the ki 
and kf pairs with high joint DOS should dominate the 
quantum interference. For energies at which the inter- 
ference fringes are prominent, the regions on CEC with 
high DOS are primarily centered around the F — K di- 
rections Therefore, three scattering wave vectors, 
labeled ql, q2 and ga, might be expected to appear in the 
power spectra. Among them, however, only q2 is along 
the f — M directions and can generate the observed stand- 
ing waves. Both qi and q^ are invisible in the power spec- 
tra. There is a simple argument to account for the disap- 
pearance of qi: the time- reversal invariance. The topo- 
logical states \k t) and \ — k I) are related by the time- 
reversal transformation: | — fc i) = T|/s f) , where T is the 
time-reversal operator. It is straightforward to show that 

i-k I \u\k ^) = -{k^\u\-k ly = -i-k I \u\k T) = 

for fermions, where C/ is a time-reversal invariant oper- 
ator and represents the impurity potential of the non- 
magnetic Ag impurity. Therefore, the backscattering 
between k and — fc is quantum-mechanically prohibited. 
Most of the observed features in the interference pat- 
tern, including the extinction of wave vector 9*3, have 
been recently well explained by a full theoretical treat- 
ment [2^ based on the T-matrix approach for multiband 
systems [l^l- In addition to the existence of standing 
waves, the absence of backscattering represents the sec- 
ond and most striking aspect of our experiment, which 
makes the topological standing waves more extraordinary 
as compared to the conventional surface states on metal 
samples [H H [ll [H • 

We can obtain the dispersion of the massless Dirac 
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FIG. 3: (a) The scattering geometry. The CEC is in the 
shape of a hexagram. The dominant scattering wave vectors 
connect two points in F— K directions on CEC. ki (pink arrow) 
and kf (red arrows) denote the wave vectors of incident and 
scattered states. q{ , 02 and ^3 (blue arrows) are three possible 
scattering wave vectors, (b) Energy dispersion as a function 
of k in the F — K direction. The data (black squares) are 
derived from FFT in Fig. 2. The red line shows a linear fit to 
the data with vp = 4.8 x 10^ m/s. The error bars represent 
the resolution of FFT (see the caption of Fig. 2). 



fermions in the F — K direction using the interference 
patterns and their Fourier transforms. For 52, the scat- 
tering geometry determines (72 = V^k [see Fig. 3(a)], 
where k is the wave vector in the F — K direction at a 
given energy. The resulting k values vary linearly with 
energy [Fig. 3(b)]. The slope of the linear fitting provides 
a measurement of the Dirac fermion velocity vp , which 
is 4.8 X 10^ m/s. In addition, the energy of the Dirac 
point is estimated to be 0.25 eV by the intercept of the 
dispersion with the energy axis. These observations are 
in agreement with the results from the first-principles cal- 
culation and the ARPES data 0, [l^] • More importantly, 
the unoccupied states, which are inaccessible to ARPES, 
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can be probed by the standing waves with STM. 

Interference fringes are also found at the step edges on 
the surface [1^ [Figs. 4(a) to 4(h)]. Similar to the case 
of Ag trimmers, the standing waves produced by steps 
are predominantly propagating along the F — M direc- 
tion. The fringes are clearly visible even at the negative 
bias voltages probably owing to the stronger scattering 
potential compared to that of the Ag trimmers. The dis- 
persion curve deduced from these patterns again shows a 
linear relation between the scattering wave vector and the 
energy [Fig. 4(i)]. Using the slope of the linear fitting to- 
gether with the same scattering geometry as that for the 
Ag trimmer, the Fermi velocity is found to be 4.8 x 10^ 
m/s, the same as that obtained from the standing waves 
caused by the Ag impurities. 




FIG. 4: Standing waves on the upper terrace by a step edge 
((a) to (h)). All the images are dl /dV maps at various bias 
voltages of an area of 35 nm by 35 nm. Imaging conditions: 
1=0.1 nA. (i) Energy dispersion deduced from the standing 
waves at the step edge. The dispersion is a function of the 
scattering wave vector q. The inserted STM image shows the 
step that produces the standing waves in (a) to (h). 

The existence of standing wave strongly supports the 
surface nature of topological states. An important issue 
that immediately arises is whether the topological states 
respond differently to the magnetic and the nonmagnetic 
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impurities. Theoretically, it was pointed out 0, |3 
that a time-reversal breaking perturbation, such as mag- 
netic impurities, can induce scattering between the states 
\k t) and — |fc i) and open up a local energy gap at the 
Dirac point. It remains an open question to observe the 
distinct signature of time- reversal breaking in topological 
insulators. 

We thank S.-C. Zhang, X.-L. Qi, Y. Ran and S.-Q. 
Shen for valuable discussions. The work is supported 



by NSFC and the National Basic Research Program of 
China. The STM topographic images were processed us- 
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Note added. At the completion of this manuscript for 
submission, we became aware of related work by P. Rou- 
san et al. [32'| . The authors reported STM study of scat- 
tering from disorder in BiSb alloy. 
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